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THEORY  OF  ELECTRON  BEAM  TRACKING 
IN  REDUCED-DENSITY  CHANNELS 


I.  INTRODUCTION 


The  propagation  of  a  charged  particle  beam  through  a  dense  gas  is 
degraded  by  collisions  with  the  gas  molecules.^  This  degradation  limits 
the  useful  range  of  charged  particle  beams,  and  thus  limits  their  utility 
in  applications  such  as  electron-beam  welding,  ’  directed -energy  weapons, 
and  inertial  confinement  fusion  reactors. One  means  of  extending  the 
range  in  a  given  environment  is  to  create  a  reduced-density  channel  using 
previous  beam  pulses  or  external  means  to  heat  the  gas.  The  heated  gas 
expands  to  reach  pressure  balance  with  the  background  gas,  leaving  a 
long-lived  rarefied  channel  that  can  dramatically  enhance  the  range  of 
subsequent  beam  pulses. 

To  benefit  from  the  channel,  the  beam  must  remain  inside.  Thus,  the 
deflection  forces  produced  by  the  channel  on  the  beam  are  important.  This 
importance  has  been  demonstrated  in  experiments  using  beams  injected  into 
hot,  highly  ionized  channels.  These  channels  quickly  expel  the  beam 
through  beam-induced  return  currents.7  Rarefied  but  highly  conducting 
channels  do  not,  therefore,  benefit  the  beam,  unless  strong  external  guide 
fields  are  provided. 

An  attractive  force  that  keeps  the  beam  in  the  channel  would  be  ideal. 

An  electrostatic  force  produced  by  low-level  channel  preionization  has  been 

q 

proposed  by  Lee,  but  this  force  is  weak,  short-lived,  partially  repulsive, 

9  10 

and  has  never  been  observed  experimentally.  ’  Much  more  important  is  the 
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magnetic  channel  force  recently  discovered  by  Welch ^  in  simulations  of 

electron  beams  propagating  through  initially  un-ionized  air.  The 

attraction  betveen  the  lov-density  channel  and  the  beam  results  from  beam 

impact  ionization  and  an  electrical  conductivity  dependent  on  the  plasma 

electron  temperature.  The  higher  temperature  in  the  rarefied  channel 

depresses  the  conductivity  and  return  current  in  the  channel,  thus  shifting 

the  centroid  of  the  net  current  toward  the  channel  axis.  Magnetic 

attraction  betveen  the  beam  current  and  net  current  pulls  the  beam  into  the 

channel.  This  attraction  has  been  observed  experimentally  by  Murphy  et 
12 

al.  using  deep,  discharge-produced  channels,  and  by  Bieniosek  using  a 
double-pulse  electron  beam,  as  reported  in  Ref.  11. 

In  this  paper  we  provide  a  theoretical  basis  for  electron  beam  guiding 
by  density  channels.  In  particular,  we  derive  analytic  formulas  for  the 
net  deflection  force  produced  by  an  initially  un-ionized,  reduced-density 
channel  on  a  rigid-rod,  azimuthally  symmetric,  ultrarelativistic  electron 
beam,  tfe  assume  for  simplicity  that  the  channel  gas  density  is 
sufficiently  high  that  a  scalar  conductivity  adequately  describes  its 
electrical  properties.  Our  analysis  shows  that  the  density-channel  force 
is  strong  and  long-lived,  and  that  it  interlaces  veil  vith  the  longitudinal 
coupling  force  binding  the  beam  body  to  the  head.  The  density-channel 
force  thus  appears  capable  of  keeping  a  (stable)  beam  in  a  channel,  without 
ancillary  fields  or  forces.  Effects  that  degrade  tracking,  including 
higher-order  chemistry  and  channel  preionization,  are  discussed  as  well. 
Numerical  simulations  using  a  nonlinear  particle  code  confirm  the  analytic 
results.  In  future  work,  we  will  present  dynamical  simulations  that 
describe  the  actual  response  of  the  beam  to  the  channel. 
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II.  GENERAL  DESCRIPTION:  THEORY  AND  EXPERIMENT 


A.  Physical  Mechanisa 

External  currents  or  charges  can  deflect  a  beam  provided  they  are 
distributed  asymmetrically  about  the  beam.  For  a  beam  passing  through  a 
neutral  gas,  the  electrical  conductivity  of  the  gas  must  be  azimuthally 
asymmetric  vith  respect  to  the  beam,  either  through  gradients  in  the  gas 
ionization  or  gradients  in  the  gas  density.  In  this  study,  ve  consider 
beam  deflection  arising  from  the  density  gradients  of  a  rarefied  channel. 
For  simplicity,  ve  assume  that  the  heavy  gas  molecules  and  ions  are 
immobile  during  each  short  beam  pulse.  Ve  additionally  assume  that  the 
dominant  source  of  ionization  is  beam  impact  ionization  (including  that 
from  the  energetic  secondaries),  and  that  it  is  a  local  and  instantaneous 
process.  Other  mechanisms  of  ionization  and  deionization  are  considered 
separately  in  Sec.  V-B. 

To  understand  hov  a  density  channel  can  deflect  a  beam,  let  us  write 
the  conductivity  of  a  weakly  ionized  gas  as 


c  * 


(1) 


where  n  is  the  number  density  of  the  plasma  electrons,  Nv^  is  their 
momentum-transfer  collision  frequency,  m  and  e  are  the  electron  mass  and 
charge,  respectively,  and  N  is  the  gas  number  density.  The  reduced 
collision  frequency,  v^,  equals  an  average  of  the  plasma-electron  speed 
times  the  cross  section  for  momentum  transfer  to  the  gas  molecules.  In  a 
veakly  ionized  gas,  is  a  function  only  of  the  plasma  electron 
temperature  Te,  and  it  rises  with  Te  for  most  gases. 

Conductivity  is  generated  as  the  beam  electrons  collide  vith  and  ionize 
the  gas  molecules.  The  collision  rate  is  proportional  to  the  gas  density 


N,  and  hence  the  degree  of  ionization,  n/N,  is  independent  of  N  and  is 
symmetric  about  the  beam.  Asymmetry  in  a  can  then  arise  only  through 
and  its  dependence  on  Tg.  Asymmetry  in  Tg  and  a  temperature-dependent 
collision  frequency  are  thus  needed  to  produce  beam  deflection. 

The  temperature  Tg  is  determined  by  the  heating  and  cooling  rates  of 

the  plasma  electrons.  The  temperature  is  high  in  low-density  regions  where 

the  collisional  cooling  rate  (proportional  to  N)  is  reduced.  In  most 

gases,  the  elevated  temperature  raises  and  lowers  a.  This  preferential 

reduction  in  plasma  conductivity  lowers  the  plasma  return  current  flowing 

in  the  channel.  The  density  depression  thus  produces  a  perturbed  channel 

current  that  flows  parallel  to  the  beam  current  and  attracts  It 

magnetically.  The  channel  force  is  attractive  and  long-lived  as  long  as 

n/N  is  symmetric  about  the  beam  and  v  rises  with  T  .  A  mathematical 

me 

formulation  of  this  model  is  given  in  Sec.  III-A.  First,  however,  we 
briefly  review  the  existing  experimental  evidence  for  the  effect. 

B.  Tracking  Experiments 

12 

Experiments  at  the  Naval  Research  Laboratory  by  Murphy  et  al.  and  at 
McDonnell-Douglas  Research  Laboratory  by  Bieniosek^*  have  confirmed  that 
electron  beams  are  attracted  to  density  channels.  Murphy  et  al.  propagated 
a  1  MeV,  10  kA,  35  ns  electron  beam  from  the  Pulserad  310  generator  into  a 
channel  formed  by  an  electrical  discharge.  The  discharge  was  directed 
through  roughly  two  meters  of  air  at  760  torr  using  a  Nd::glass  laser  to 
prepare  a  weakly  preionized  path.  The  discharge  could  be  displaced  up  to 
three  cm  from  the  beam  axis,  and  it  produced  a  channel  that  grew  to  several 
cm  in  radius  at  pressure  equilibrium.  The  effect  of  the  channel  on  the 
electron  beam  was  ejection  from  hot  channels  and  attraction  toward  cool 
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channels.  This  behavior  is  expected  because  hot  channels  are  highly 

conducting  vhile  cool  channels  are  not.  Although  the  tracking  effect  vas 

unambiguous,  precise  quantitative  measurements  were  hampered  by  large  hose 

growth  and  uncertainties  in  the  channel  parameters.  Continuing  experiments 

13 

are  being  performed  using  the  5  MeV  SuperlBEX  generator. 

The  experiments  of  Bieniosek**  were  performed  on  the  two-pulse  MEDEA  II 
device.  The  parameters  of  each  pulse  were  similar  to  those  of  Pulserad, 
and  a  variety  of  pulse  separation  times,  gas  mixtures  and  gas  densities 
were  employed.  The  first  pulse  was  deflected  by  a  transverse  magnetic 
field  that  was  turned  off  prior  to  injection  of  the  second  pulse. 

Attraction  of  the  second  pulse  to  the  density  channel  formed  by  the  first 
pulse  demonstrated  the  effect.  Since  the  density  reduction  was  modest,  the 
attraction  was  much  weaker  than  in  the  Pulserad  experiments. 


III.  TRACKING  FORMALISM 


A.  Force  Equation 

In  this  section  we  derive  formulas  for  the  magnetic  deflection  force 
produced  by  an  initially  un-ionized  density  channel  on  an  ultrarelativistic 
electron  beam.  The  force  is  termed  tracking  if  it  pulls  the  beam  toward 
the  channel,  and  detracking  if  it  pushes  the  beam  away.  To  simplify  the 
analysis,  we  assume  that  the  beam  is  azimuthally  symmetric  about  a  common 
axis,  and  we  Ignore  axial  and  temporal  variations  relative  to  transverse 
variations  to  order  3/3z,  3/3C  «  l/r^;  here  r^  is  a  characteristic  beam 
radius,  and  C  *  ct-z  measures  distance  behind  the  beam  head.  The  force 
equation  then  reduces  to  a  static,  two-dimensional  calculation,  with  the 
beam  and  plasma  treated  as  infinitely  long  at  each  location  £.  To  this 
order  of  approximation,  there  is  no  transverse  current  flow.  Furthermore, 
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because  an  infinitely  long,  azimuthally  symmetric  ring  of  beam  current 
produces  no  fields  or  forces  inside  itself,  the  ring  can  feel,  by 
reciprocity,  no  net  force  from  currents  and  charges  inside.  Hence,  only 
currents  and  charges  outside  the  ring  can  deflect  it.  See  Fig.  1. 

Consider  now  a  plasma  return  current  of  density  oE  ,  where  E  is  the 
axial  electric  field.  In  the  direction  linking  the  beam  and  channel  axes, 
0  *  0,  this  current  magnetically  deflects  a  beam  ring  of  radius  r  and 
current  5lb  with  a  distributed  net  force  given  by  the  Biot-Savart  law: 


F«<r.O  dr 


Slb(r,0  Jdr' 
r 


r' 


2n 

d0  cos0 

% 


(2) 


where  8Ib(r,C)  *  2nr  dr  jb(r,C)  and  jb  is  the  beam  current  density. 
Summing  over  all  beam  rings  and  integrating  by  parts  yields  the  average 
force  per  beam  electron: 


Ft»<«  - 


r  F„<r-« 


•  3i.  *  2" 
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r  0 
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where  Ib  is  the  total  beam  current  and 


ib(r,0 


r 

Jdr  2nr 
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jb(r.O 

Ib(C) 


(3) 


(4) 
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is  the  fraction  of  beam  current  within  radius  r.  In  deriving  Eq.  (3),  we 
have  invoked  definition  (1)  for  a  and  assumed  that  n/N  is  symmetric  about 
the  beam. 

As  shown  in  Sec.  IV-A,  the  field  E  soon  becomes  azimuthally  symmetric, 

Z 

and  hence  it  can  be  pulled  outside  the  integral  over  0.  We  can  then 
rewrite  the  angular  integral  by  invoking  the  assumption  that  \>m  depends  on 
6  only  through  its  dependence  on  N: 


jde  ESS  .  .  jde  slne  |g  (M 


=  Jde  sine  ln(N0/N), 


where  N  is  the  ambient  gas  density  outside  the  channel,  and  we  define 


3(v  /\>  ) 
v  mo  nr 

31n(NQ/N)  ‘ 


Here  v  is  a  characteristic  value  of  the  reduced  collision  frequency, 
mo 

we  now  assume,  as  discussed  in  Sec.  V-A,  that  q  is  a  gas  parameter 
independent  of  (r,6),  we  can  rewrite  force  equation  (3)  as 


Ftm<«  ’  Idr  *b  °oZz  JdS  sl"0  k  ln(No/N>  ' 


where  oQ(r,0  *  e  n'/ra^vmo' 


B.  Linearized  Analysis 

Consider  first  a  small  channel  offset,  yc,  satisfying 
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yc  <<  V  V 


(8) 


where  r.  is  the  beam  radius  and  r  is  the  channel  radius.  We  can  then 
D  c 

expand  the  gas  density  about  the  beam  axis  as 

ln[N0/N(r,e)]  *  In  [NQ/N(r)]  -  yccose  |p  In  [NQ/N(r)]  +  . . .  (9) 

This  approximation  reduces  force  equation  (7)  to 


Ftm<C>  -  2Jr  Wc  Jdr  ”oEz  I?  1”<VN>  '  <10> 

0 

with  the  solution  determined  by  the  radial  distributions  of  the  beam 

current,  the  plasma  current,  and  the  density  channel. 

A  particularly  simple  example  is  a  square  density  channel  for  which 

N  >  N  for  r  <  r  and  N  *=  N  for  r  >  r  .  Then, 
c  c  o  c 

Ft.(C)  »  ^,yc  ln(H0/Nc)  f(VC>  jb(rc,0  i„<V0  ,  (11a) 


where  the  plasma  return-current  fraction,  f  »  -a  E  /jK,  is  positive  for 

O  Z  D 

rising  or  constant  If  the  beam  and  its  ionization  take  self-similar 

2  2  2  2 

Bennett  profiles  of  the  form  j^(r,C)  *  1^(0^  /n(r  +r^  )  and  i^(r,C)  = 
r2/(r2+rb2) ,  Eq.  (11a)  reduces  to 


2eIb  yc 

Ft*<«  '  ,£  V  ^  ln(VNc> 


r2/r„2 


(!  +  rc2/rb2)3 


(lib) 


Several  features  can  be  deduced  from  these  results.  First,  the 

tracking  force  can  be  large  and  even  exceed  the  magnetic  self-pinch  force, 

F  *  e(l-f )Iu/r.  c,  that  holds  the  beam  together.  Although  the  beam  is 
o  0  0 
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likely  to  distort  azimuthally  when  Ftm  >  Fq,  such  distortion  should  not 

disrupt  channel  tracking.  Moreover,  Ftffl  «  Fq  in  most  applications,  and 

the  rigid-rod  treatment  is  then  well  justified.  Second,  the  sign  of  F^  is 

tm 

independent  of  the  relative  size:  of  the  beam  and  channel,  and  it  is 

tracking,  Ftm/yc  >  0»  provided  q,f  >  0  and  Nq  >  N£.  Third,  F  has  a  broad 

maximum  when  the  beam  and  channel  radii  are  comparable,  but  diminishes 

rapidly  as  the  channel  becomes  much  smaller  or  larger  than  the  beam. 

Fourth,  F  is  long-lived,  persisting  as  long  as  plasma  return  current 

flows  at  the  channel  edge.  And  fifth,  F  depends  logarithmically  on  th« 

channel  gas  density.  This  weak  dependence  suggests  that  deep  channels  are 

not  needed  to  effect  strong  tracking.  Only  in  shallow  channels,  (N  -  N  ) 

«  N  ,  is  the  force  is  sensitive  to  the  channel  depth,  F..  *  (N  -  N  ). 

o  tm  o  c 

For  small  offsets,  the  tracking  force  varies  linearly  with  y  .  This 
suggests  that  the  characteristic  axial  distance  needed  to  pull  the  beam 
into  the  channel  is  given  by  the  harmonic-oscillator  result: 


z 


t 


n 

2 


ymc 


F  /y 
.  tm  •'cj 


1/2 


(12) 


where  y  is  the  relativistic  mass  factor  for  the  beam.  Using  force  equation 

(lib), 


z 


t 


(1  ♦  rc2/rb2)3 


2q  ln(No/Nc) 


1/2 


(13) 


1/2 

where  =  2nrb[IA/(l-f >1^]  is  an  average  betatron  oscillation 

1  3 

wavelength  for  the  beam  electrons,  and  1^  =  ymc  /e  =  17y  kA  is  the  beam 

Alfven  current. 
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C.  Nonlinear  Analysis 

A  linear  analysis  simplifies  the  mathematics  but  is  restricted  to  small 
offsets  yc-  For  completeness,  we  present  a  nonlinear  analysis  valid  for 
arbitrary  offset.  Ve  also  choose  a  more  realistic,  rounded  profile  for  the 
channel  gas  density  that  is  amenable  to  numerical  treatment: 


N(r')  =  N 


2  1p 
(r'/rcJ  +  6 

(r'/rc)2  +  1 


2  2  2 

where  r'  *  r  +  y  -  2ry  cos9,  and  &  <  1  and  p  >  0  are  constants.  In  the 

c  c 

channel  center,  N(0)  =  N£  =  No8p. 

Profile  (14)  is  easily  differentiated  in  force  equation  (7): 


a  a  [r2  +  yc2-  2ycr  cos®  +  r  2 

§9  in(N  /N)  =  pfem-j - V“T - : - 

r  +  y„  -  2y  r  cos0  +  rS 
c  c  c 


2pycr  sin© 

**2  ?  2 
r  +  y„  -  2y  r  cos©  +  r 
c  c  c 


2py  r  sin© 
_  c _ 

2  2 

r  +  y  -  2y  r  cos© 
c  c 


2  ’ 
rc  6 


The  integral  over  ©  can  then  be  performed  using  the  general  result  that 


,  2 
sin  6 

a  -  b  cos© 


a  -  Ja2  -  b2 


for  a  >  b.  After  considerable  algebra,  this  reduces  Eq.  (7)  to 


p,*<«  -  -  I  T  1 1  Vs,  l(1  -  S)rc2 
c  0  v 

♦  [r4  *  2r2(r.2J  -  yc2)  ♦  <yc2.  rc2S)2]1/2 

-  fr*  *  2'V-  »,’>  *  <»«*♦  r=2)2]1/2)  ■ 
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Consider  now  a  beam  vith  a  Bennett  profile  of  constant  radius  r^.  The 
degree  of  .onization,  n/N,  produced  by  impact  ionization  is  then 
self-similar  vith  the  beam.  Ve  also  assume  that  Ez  is  independent  of 
radius  r,  which  is  usually  a  fair  approximation  as  discussed  in  Sec.  IV-A. 
The  assumption  that  q  and  are  independent  of  (r,6)  implies  that  the 
reduced  collision  frequency  takes  the  peculiar  form 


vm(r,0,O  = 


v(C) 

mo 


-  ln(lNo/N(r,e>]q)  ‘ 


This  form  is  valid  only  if  q  <  1  and  the  minimum  gas  density  exceeds 
Nq  exp(-l/q).  The  validity  of  Eq.  (18)  and  typical  values  for  q  and 
are  discussed  in  Sec.  V-A. 

Vith  these  assumptions,  Eq.  (17)  can  be  integrated.  The  result  is 


F,m<«  -  8?  h  k>  -  *<s>]  - 


where  Ip  is  the  total  plasma  return  current  and  k  and  g  are  defined  through 
the  relationships 


(2r  r  y  )2  y  4-  (r  2-  r  V 

g(S)  «  - -  -ft-  S  ln[a(S)]  +  2 - 

d  (8)b  d(8)4 


(20a) 


d(S)  =  [(yc2+  rc28)2  +  2rc2(yc2-  rc2S)  +  rb4]  , 


(20b) 


2  „  2 


|d(t)2.  rc28  -  yc2-  rb2)rb2 

rw, 2  _  2 ..2  .  2 


(y„  +  r /«)d(S)%  (y_  +  r  6)  *  (y_  -  r/8)r 


2  2  ’  (20c) 


k(8)  «  1  +  P  [ln(8)  +  gx(l)  -  gl(8)]  , 


(20d) 
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and 


y  +  r  S  +  r, 

gi(S)  =  — - S - “  ln[«(6)].  (20e) 

1  d(6)2 

Equation  (20d)  for  the  normalization  coefficient  k  vas  obtained  using  the 
definition  of  the  plasma  return  current, 


(21) 


where  n/N  has  the  same  Bennett  profile  as  the  beam,  and  E  is  again  assumed 

Z 

to  be  independent  of  r. 

The  nonlinear  result  (19)-(20)  offers  generality  at  the  expense  of 
complexity.  The  principal  additional  information  is  that  F  peaks  when 
the  beam  lies  at  the  channel  edge  (where  31n[No/Nj/3r'  peaks),  and  that  F 
remains  tracking  and  falls  off  nearly  quadratically  with  yc  when  y£  »  rc« 
(The  fall-off  is  faster,  however,  if  the  beam  or  channel  is  sharp-edged; 
for  example,  Ftm  equals  zero  if  the  beam  is  entirely  outside  the  channel.) 
The  nonlinear  result  also  confirms  that  Ftm  rises  nearly  linearly  with  the 
chemistry  parameter  q  and  weakly  with  the  channel  depth.  These 
dependencies  are  shown  in  Fig.  2. 


IV.  FIELD  EQUATIONS 

A.  Neglecting  the  Dipole  Electric  Field 

The  preceding  tracking  analysis  is  predicated  on  the  assumption  that 
the  electric  field  is  symmetric  about  the  beam  and  varies  weakly  with 
radius  r.  Here  we  justify  that  assumption  in  the  high-conductivity  region 
where  the  plasma  return  current  and  magnetic  tracking  are  important.  In 
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this  region,  the  electrostatic  potential  and  transverse  electric  field  are 
shorted  out,  leaving  an  axial  electric  field  given  by 

Eg  -  -  I?  •  (22) 

The  axial  vector  potential  A  and  field  E  are  determined  by  the  Lorentz 

z 
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wave  equation,  given  in  the  frozen-field  limit  by 

7!2a  -  -  T  [jb  ♦  <**]  •  <23> 

The  behavior  of  E  is  most  easily  analyzed  using  the  linearization 

Z 

employed  in  Sec.  III-B.  For  small  channel  offsets,  the  vector  potential 
can  be  expanded  in  terms  of  monopole  and  dipole  components: 


A(r,0,t)  =  AQ(r,C)  +  Aj(r,0  cos0  +  ... 


(24) 


We  again  consider  ohmic  heating  only,  but  for  convenience,  ve  replace  Eq. 
(18)  with  the  modified  form, 


vm(r,e,t) 


v 

mo 


T5z(r,e,t)  N 

N(r,9)  E 


(25) 


where  q  is  constant  and  Eq/No  is  a  reference  field  strength.  Imposing 
expansion  (9)  for  N  and  expansion  (24)  for  A  on  Eqs.  (22),  (23)  and  (25) 
produces,  after  some  manipulation,  the  linearized  equations 


and 


» 


3_  1  3_ 
3r  r  3r 


r 


-!fo 

3C 


f?  ln<VN>  » 


(26) 


(27) 
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vhere  the  monopole  conductivity  a( r,C)  is  computed  using  E  *  -3A  /3C  in 

z  o 

Eq.  (25)  for  vffl. 

The  characteristic  relaxation  length  for  the  monopole  potential  Aq  can 
be  determined  from  Eq.  (26).  For  self-similar  Bennett  profiles,  this 
length  is  given  by 


no(0,Orb' 


In  1  +  b  /r. 


(28) 


where  o(0,Q  is  the  on-axis  monopole  conductivity,  and  b  »  r^  is  the 
radius  within  which  space-charge  neutralization  occurs.  The  monopole 
potential  in  this  case  is  given  by 


Ao(r,C) 


(29) 


where  f  is  the  plasma  return-current  fraction.  The  weak,  logarithmic 
dependence  on  r  justifies  treating  the  monopole  electric  field,  -3Aq/3C,  as 
independent  of  r  in  the  tracking  analysis. 

Similarly,  the  relaxation  length  of  the  dipole  potential  A^  can  be 
determined  from  Eq.  (27),  and  for  self-similar  Bennett  profiles  is  given  by 

_  no(0,Orb2 

<ml  -  <1^)  - 2c  ~  *  (30) 

This  equation  was  first  derived  by  Slinker  et  al.*3  in  studies  of  hose 
instability.  According  to  Eqs.  (28)  and  (30),  the  dipole  decay  length 
is  roughly  an  order  of  magnitude  smaller  than  the  monopole  decay  length 
C|nQ.  Neglecting  the  dipole  electric  field,  -3A^/3C,  relative  to  the 
monopole  field,  -3Aq/3C»  is  thus  well  justified  in  the  tracking  analysis. 

Ve  state  without  proof  that  the  dipole  field  acts  to  reduce  the  channel 
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force,  and  that  the  linearized  tracking  result  (10)  can  be  recovered  from 
Eq.  (27)  by  dropping  9A^/dC  and  using  the  fact  that  the  local  magnetic 
deflection  force  equals  -eVjA. 


B.  Longitudinal  Coupling  and  the  Resistive  Hose  Instability 

The  channel  tracking  force  calculated  in  the  previous  sections  varies 
vith  distance  t  into  the  pulse.  As  the  beam  bends  in  response  to  this 
force,  it  perturbs  the  electromagnetic  fields.  The  plasma  retards  the 
growth  and  decay  of  these  perturbations  and  thus  generates  a  longitudinal 
coupling  force,  as  the  fields  produced  by  early  beam  segments  persist  to 
guide  later  segments.  Although  a  detailed  analysis  of  dynamical  effects  is 
beyond  the  scope  of  this  paper,  the  importance  of  coupling,  both  to 
tracking  and  to  the  resistive  hose  stability,  varrants  a  brief  discussion. 
Further  discussion  of  coupling  is  given  in  the  Appendix. 

Coupling  develops  from  dipole  electric  fields  that  are  short-lived 


«w  «  Cmo)  relative  to  the  monopole  fields  responsible  for  channel 
tracking.  As  a  consequence,  the  channel  force  exceeds  the  coupling  force 


as  long  as  the  channel  offset  and  size  remain  modest,  y£  <  r£  -  r^,  and  as 
long  as  the  simple  chemistry  model  discussed  earlier  remains  valid.  As 


discussed  in  Sec.  V-A,  however,  the  latter  assumption  fails  at  late  C,  as 


higher-order  chemistry  causes  the  channel  force  to  change  sign  and  repel 
the  beam.  The  beam  tail  would  then  be  ejected  from  the  channel,  even 
though  the  head  tracks,  in  the  absence  of  coupling. 

Although  coupling  is  relatively  unimportant  early  in  the  beam  head,  it 
becomes  important  in  the  body  where  cumulative  ionization  produces  a  large 


dipole  decay  length  If  the  channel  force  changes  sign  quickly  over 


*ml’ 


the  coupling  force  will  dominate  and  cause  subsequent  portions  of  the 
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beam  to  follov  the  head,  regardless  of  the  channel  force.  Coupling  can 

thus  keep  the  entire  beam  in  the  channel,  provided  the  channel  force  is 

tracking  initially  and  does  not  become  detracking  until  is  large. 

Although  coupling  keeps  the  beam  together,  it  also  drives  the  resistive 

hose  instability,  as  field  perturbations  become  out  of  phase  with  beam 
1  /  1 

oscillations.  *  Because  the  phase  lag  and  growth  rate  of  this 
instability  are  governed  by  the  magnetic  dipole  decay  length  slow 
changes  in  the  channel  force  over  should  not  excite  strong  hose  growth. 
To  the  contrary,  a  channel  force  that  is  tracking  and  stronger  than  the 
coupling  force  should  reduce  hose  growth  by  providing  a  fixed  central  axis 
of  attraction.  Because  hose  is  a  convective  instability  and  because 
rises  with  t,  damping  is  most  important  in  the  beam  head1^  which  is  where 
density  tracking  is  strongest.  Density  tracking  can  thus  suppress  the  hose 
instability,  as  well  as  keep  the  beam  in  the  channel. 

C.  Electrostatic  Forces 

The  preceding  magnetic  tracking  analysis  is  restricted  to  the  high- 
conductivity  region  where  the  electrostatic  fields  and  forces  are  small. 

The  electrostatic  fields  are  not  small,  however,  early  in  the  beam  head 
where  the  conductivity  is  low.  Nevertheless,  several  effects  conspire  to 
reduce  the  impact  of  these  fields  on  beam  steering. 

Electrostatic  deflection  forces  arise  as  the  channel  conductivity 
asymmetrically  neutralizes  the  electrostatic  field  of  the  beam,  thereby 
polarizing  the  plasma  charge.  Because  deflection  develops  from  charge 
polarization  rather  than  from  a  unipolar  charge,  the  deflection  force 
varies  with  radius  r,  and  hence  its  average  value  is  small. 
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A  second  mitigating  effect  is  that  the  plasma  charge  polarizes  on  a 
monopole  decay  length  given  by 

Ceo  *  c/4nff*  (31a) 

but  relaxes  on  a  dipole  decay  length  given  by 

tel  =  c/(l-q)2no.  (31b) 

These  relationships  can  be  derived  using  charge  conservation  and  Gauss's 

o 

lav.  Because  the  two  lengths  are  comparable,  electrostatic  deflection 
decays  nearly  as  rapidly  as  it  arises,  thus  reducing  its  peak  value.  By 
contrast,  magnetic  deflection  attains  full  value  because  it  decays  slowly 
relative  to  its  growth, 

A  third  limitation  is  the  rapid  growth  of  o  from  beam  ionization.  This 
ionization  quickly  reduces  *tel»  thereby  confining  electrostatic  deflection 
to  near  the  beam  head.  Because  longitudinal  coupling  is  not  well  developed 
there,  electrostatic  deflection  is  ineffectual  at  steering  the  beam,  as 
well  as  weak. 


V.  GAS  CHEMISTRY 


A.  Plasma-Electron  Collision  Frequency 

The  chemistry  parameter  q  (or  q)  determines  the  sign  and  magnitude  of 
the  channel  force,  and  it  thus  largely  determines  the  behavior  of  beams  in 
density  channels  in  a  given  gas.  In  this  section,  we  discuss  the  physical 
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basis  for  q  and  vmo,  and  determine  their  dependence  on  gas  chemistry  and 
beam  heating. 

As  noted  earlier,  the  reduced  collision  frequency  in  a  weakly 

ionized  gas  equals  an  average  of  the  plasma  electron  speed  times  the  cross 

section  for  momentum  transfer  to  the  gas  molecules.  The  average  electron 

speed  increases  as  the  square  root  of  the  plasma  electron  temperature  Te, 

and  for  most  gases  this  dependence  causes  v  to  rise  with  T  .  If  v  rises 

m  e  m 

with  T  ,  and  if  T  falls  with  gas  density  N,  the  parameter  q  in  Eq.  (6)  is 

c  c 

positive,  and  the  channel  force  is  tracking.  Otherwise,  the  force  is 
de tracking. 

The  dependence  of  Tfi  on  N  can  be  determined  as  follows.  The  plasma 
electrons  gain  energy  from  both  ohmic  heating  and  direct  beam  heating,  and 
lose  energy  via  collisions  with  the  gas  molecules.  Ohmic  heating  derives 
from  the  electric  fields  of  the  beam  and  plasma,  while  direct  beam  heating 
is  a  by-product  of  beam  impact  ionization.  Setting  the  cooling  rate  of  the 
plasma  electrons  equal  to  their  heating  rate  yields 

nNvc  .  cE2  +  fh  (jb/e)  N  jjf  ,  (32a) 

where  is  the  reduced  inelastic  cooling  rate  for  the  plasma  electrons, 
de/dx  is  the  reduced  energy-loss  function  for  the  beam,  and  f^  is  the 
fraction  of  deposited  beam  energy  going  into  heating  of  the  plasma 
electrons.  For  a  given  gas,  is  a  function  of  Tfi  only,  and  de/dx  is  a 
function  of  beam  energy  r  only.  For  relativistic  electron  beams  in  air, 
de/dx  -  2-3  keV/cm-atm,  and  f^  ~  0.2. 

Equation  (32a)  can  be  rewritten  by  dividing  by  nN  and  substituting 
Eq.  (1)  for  cr.  This  yields  a  transcendental  equation  for  the  plasma 
electron  temperature  Te: 
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1 

N  * 


(32b) 


v  (T  ) 
c'  e' 


“W 


f±  de 
e  dx 


n/N 


\  i 


To  lowest  order,  the  electric  field  E  and  the  degree  of  ionization  n/N  are 
independent  of  gas  density  N.  As  a  consequence,  lowering  N  raises  the 
heating  rate  on  the  right-hand  side  of  Eq.  (32b).  Tg  is  thus  elevated  in 
regions  of  low  N.  The  elevation  is  usually  modest,  however,  because  vc 
rises  rapidly  with  Tfi  for  temperatures  of  interest,  Tg  -  1  eV.  The  weak 
dependence  of  Tg  on  N  is  responsible  for  the  weak  dependence  of  the 
tracking  force  on  N. 

The  chemistry  parameter  q  is  determined  by  Eqs.  (6)  and  (32b)  and  the 

dependence  of  v  and  v  on  T  .  A  much  simpler  estimate  can  be  obtained, 
m  c  e 

however,  by  ignoring  direct  beam  heating  relative  to  ohmic  heating.  The 

reduced  collision  frequency  \»  can  then  be  expressed  as  a  function  of  E/N, 

without  reference  to  Tfi.  The  neglect  of  direct  beam  heating  is  usually 

valid  when  the  beam  and  its  fields  are  intense. 

As  an  example,  we  plot  in  Fig.  3  the  reciprocal  of  vm  as  a  function  of 

the  reduced  field  E/N  in  air.  This  curve  was  computed  from  the  electron 

drift  speed,  w  =  eE/mN\>m,  as  calculated  by  Heylen*7  and  tabulated  by 
18 

Dutton.  For  convenience,  E/N  is  normalized  to  the  avalanche  air 

breakdown  strength,  E^/N  *  25  kV/cra-atm.  Much  above  this  field,  avalanche- 

induced  detracking  is  likely,  as  discussed  in  the  next  section. 

Figure  3  indicates  that  the  chemistry  model  (18)  is  a  crude  but 

adequate  approximation  over  limited  ranges  of  E/N.  In  this  model,  the 

parameters  q  and  vmo  are  weak  functions  of  the  electric  field  E(C),  which 

is  assumed  to  be  independent  of  r.  Typical  values  for  air  are  q  s  0.2  and 

v  »  10'7  cm^/s  for  0.1  E.  /N  <  E/N  <  E. /N;  at  lower  fields,  q  «  0.5. 
mo  d  ~  D 
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Replotting  to  fit  form  (25)  shows  equally  good  agreement,  with  q  a  0.2 
at  high  fields  and  q  =0.5  at  the  lover  fields.  All  of  these  parameters 
are  insensitive  to  the  nature  of  the  heating  mechanism,  provided  the 
cooling  rate  vc  remains  a  strong  function  of  T  .  The  addition  of  direct 
beam  heating  should  not,  therefore,  notably  affect  q,  q,  or  v  . 

B.  Other  Sources  and  Sinks  of  Ionization 

Density  tracking  derives  from  simple  but  subtle  chemistry  properties, 

and  this  raises  the  concern  that  other  processes  may  be  equally  important. 

In  this  section,  we  examine  the  effects  of  higher-order  chemistry  on  n/N. 

19 

A  more  careful  numerical  analysis  was  performed  by  Keeley  using  a 
detailed  air  chemistry  code. 

The  behavior  of  n/N  can  be  examined  using  a  generalized  continuity 
equation  for  the  plasma  electrons: 

||  -  Ns.jb/e  -  <va  -  0a)nN  -  0rn2  -  W-(nv)  ,  (33) 

where  s.  is  an  effective  beam  ionization  cross  section,  v  is  an  electron 

J.  21 

avalanche  coefficient,  6  is  an  electron  attachment  coefficient,  8  is  an 

©  I 

electron-ion  recombination  coefficient,  and  w  is  the  plasma-electron  fluid 
velocity.  In  the  preceding  analysis,  all  terms  on  the  right-hand  side  of 
Eq.  (33)  were  neglected.  We  separately  assess  each  of  these  terms,  using 
the  frozen  approximation  that  3/3t  =  c  3/3C. 

The  electron  avalanche  coefficient  v  is  a  strong  function  of  E/N,  and 
thus  it  preferentially  produces  high  n/N  and  high  a  in  a  low-density 
channel.  High  channel  conductivity  increases  the  plasma  return  current  in 
the  channel,  which  repels  the  beam.  Avalanching  must  therefore  be  kept 
small  relative  to  beam  ionization  to  keep  the  net  channel  force  tracking.11 
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This  requires  high  channel  gas  density  N£  and  modest  electric  field  E.  A 
generally  sufficient  condition  is 

E/N,  <  3(Eb/N),  (34) 

where  E^/N  is  the  minimum  reduced  field  needed  to  initiate  avalanching. 
Condition  (34)  bounds  the  minimum  channel  depth  for  a  given  beam. 

The  effect  of  electron  attachment,  which  converts  free  electrons  to 
essentially  immobile  negative  ions,  depends  on  the  dependence  of  the 
attachment  coefficient  0a  on  temperature  Tg  and  gas  density  N.  If  there  is 
little  or  no  temperature  dependence,  the  effect  is  detracking,  because 
attachment  lowers  n/N  more  outside  the  density  channel  than  inside,  so  that 
plasma  return  current  flows  preferentially  in  the  channel.  This  detracking 
becomes  important  at  distances 

t  >  c/f$aNo  (35) 

from  the  head. 

Electron-ion  recombination  is  detracking  for  the  same  reason.  Namely, 
it  lowers  n/N  more  outside  the  channel  than  in.  Recombination  becomes 
important  at 

t  >  He/Vi^h  ,  (36) 

where  Nj^  is  to  be  evaluated  at  its  maximum  value. 
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C.  Plasma  Dynamic  Effects 

In  addition  to  the  chemistry  issues,  the  analysis  neglects  several 
plasma  dynamic  effects.  One  of  these  is  plasma  electron  convection 
represented  by  the  last  term  in  Eq.  (33).  Convection  perturbs  n/N  much 
like  a  source  or  sink  of  ionization,  but  can  usually  be  neglected  because 
the  electron  fluid  velocity  is  typically  small, 

v/c  <  10"3.  (37) 

A  related  dynamical  effect  is  electron  pressure  in  Ohm's  law.  Electron 
pressure  can  be  neglected  provided 

Te  «  eErb,  (38) 

which  is  usually  well  satisfied  except  for  very  narrow  beams. 

A  third  dynamical  effect  is  the  finite  time  taken  for  plasma  heating. 
This  time  can  be  neglected  provided  the  beam  parameters  vary  slowly.  A 
sufficient  condition,  based  on  ohmic  heating  alone,  is 

c  §£  «  eEw/Tg.  (39) 

The  heating  time  for  intense  beams  in  air  is  typically  sub-nanosecond. 

A  fourth  effect  is  diffusion  from  beam  ionization.  Beam  ionization 

liberates  energetic  secondaries  that  produce  additional  ionization  through 

a  cascade  process.  The  total  ionization  is  thus  spread  in  space  and  time 

20  21 

inversely  proportional  to  N.  ’  This  spread  preferentially  reduces  the 
formation  rate  of  a  in  the  low-density  channel,  an  effect  that  should 
enhance  tracking. 
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A  fifth  neglected  effect  is  momentum  transfer  between  the  plasma 

electrons  and  ions.  These  Spitzer  collisions  become  important  at  high 

_2 

degrees  of  ionization,  n/N  >  10  ,  and  cause  vm  to  decrease  rather  than 

increase  with  Tg,  thereby  producing  detracking.  In  initially  un-ionized 
gases,  however,  n/N  is  initially  small,  and  hence  collisions  with  ions  are 
usually  negligible  in  the  beam  head  where  tracking  is  most  important. 

D.  Channel  Preionization 

The  assumption  that  the  channel  is  un-ionized  prior  to  beam  injection 
fails  if  the  channel  is  deep  and  in  thermodynamic  equilibrium.  In 
equilibrium,  a  deep  channel  must  be  hot  to  maintain  pressure  balance  with 
the  background  gas,  and  hence  it  must  be  thermally  ionized.  Return 
currents  induced  in  this  ionization  can  eject  the  beam.^  In  ‘his  section 
we  determine  how  much  preionization  can  be  tolerated,  and  the  limitations 
that  thermal  preionization  imposes  on  channel  depth. 

Channel  preionization  produces  magnetic  and  electrostatic  deflection 
forces.  As  discussed  earlier,  the  magnetic  forces  develop  on  the  dipole 
decay  length  and  relax  on  the  monopole  decay  length  CmQ,  while  the 
electrostatic  forces  develop  on  the  monopole  decay  length  Ceo  and  relax  on 
the  dipole  decay  length  All  the  forces  relax  more  slowly  than  they 

develop. 

At  low  levels  of  preionization,  initially  exceeds  C|n0,  and  thus  the 
electrostatic  forces  initially  dominate.  These  forces  are  relatively 
ineffectual,  however,  at  steering  the  beam  as  discussed  in  Sec.  IV-C. 

By  contrast,  at  high  levels  of  preionization,  everywhere  exceeds 
and  the  magnetic  forces  then  dominate  everywhere.  According  to  Eqs. 
(30)  and  (31),  exceeds  **el  for  all  C  provided 
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(40) 


(l-tOna^/c  >  1, 


where 
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2 
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is  the  conductivity  resulting  from  the  channel  preionization,  n^.  Yu 

derived  a  variation  of  condition  (40)  and  termed  it  the  channel  overheat 

condition  to  signify  the  onset  of  strong  magnetic  detracking. 

In  a  density  channel,  account  must  be  taken  of  the  density  force  F  . 

To  avoid  repetition,  we  will  not  compute  the  detracking  force  from  channel 

preionization,  but  instead  determine  only  when  it  is  cancelled  by  F  . 

For  simplicity,  consider  a  square  density  channel  of  radius  r^  that  is 

uniformly  preionized  and  entirely  inside  a  square  beam  of  radius  r^  >  r£. 

Assuming  beam  impact  ionization  only,  the  uniform  channel  preionization, 

nj/N,  is  supplemented  by  uniform  beam  ionization,  n(C)/N.  The  gas-density 

depression  then  perturbs  the  plasma  current  by  an  amount  proportional  to 

the  change  in  v  , 
m 


d 


®  k<N>  "  vU  ’ 


(42) 


while  the  channel  preionization  perturbs  the  plasma  current  by 


v  (N) 


(43) 


The  current  produces  tracking  while  produces  detracking. 
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Because  the  tvo  currents  have  the  same  spatial  distribution,  the  net 
magnetic  channel  force  equals  zero  when  they  cancel,  Sj^+  Sj^  *  0.  The 
transition  from  detracking  to  tracking  thus  occurs  at  the  point  Ct  where 


n(Ct) 


ni 


V(N  ) 
m  o 


VN>  -  vv 


1  -  q  ln(N0/N) 
*  ni  q  ln(NQ/N) 


(44a) 


(44b) 


Form  (44a)  is  based  on  chemistry  relationship  (18)  for  v^,  while  form  (44b) 
is  based  on  chemistry  relationship  (25).  For  deep  channels,  the  force 
typically  transitions  from  detracking  to  tracking  once  n(t)  >  nj. 

As  discussed  earlier,  the  magnetic  coupling  force  relaxes  on  the  dipole 
decay  length  Coupling  is  thus  strong  at  Ct  provided 


(l-q)ita.rb2/c  >  C.%. 


(45) 


This  condition  is  stronger  than  the  channel  overheat  condition  (40), 
provided  Ct  >  r^.  The  latter  proviso  is  in  fact  fundamental  to  the 
analysis  which  presumes  that  the  beam  and  plasma  vary  slowly,  r^  3/9C  «  1. 
The  proviso  fails  only  for  extremely  intense  beams,  for  which  displacement 
currents  and  like  effects  must  be  incorporated. 

For  beams  of  interest,  condition  (45)  is  sufficient  and  generally 
necessary  for  channel  preionization  to  eject  the  entire  beam.  The  head  is 
ejected  because  the  net  channel  force  is  detracking,  while  the  body  is 
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ejected  because  it  is  coupled  to  the  head.  I£  condition  (45)  is  not 
satisfied,  the  detracking  portion  of  the  head  is  still  ejected,  but  later 
portions  of  the  beam  can  track  the  channel.  As  discussed  in  the  Appendix, 
detracking  condition  (45)  is  imprecise  because  it  does  not  account  for  the 
different  parametric  dependencies  of  the  channel  and  coupling  forces.  The 
Appendix  also  shows  that  condition  (45)  can  be  difficult  to  meet. 

Although  conditions  (40)  and  (45)  differ  significantly  regarding  n^/N, 

the  difference  is  often  insignificant  regarding  allowed  channel  depth. 

This  is  because  n^/N  is  a  strong  function  of  N  for  a  thermalized  but  weakly 

ionized  gas  in  pressure  balance.  To  show  this,  consider  the  Saha  equation 

which  states  that  the  degree  of  thermal  preionization  n^/N  in  a  weakly 

ionized  gas  is  proportional  to  exp(-V^/2T  ),  where  is  the  gas  ionization 

potential  and  T  is  the  local  gas  temperature.  The  gas  is  weakly  ionized 
© 

only  if  Tg  «  W^/2,  in  which  case  nj/N  is  a  very  strong  function  of  Tg.  In 
pressure  balance,  the  gas  density  varies  inversely  with  temperature, 

N  «  T  so  that  ni/N  varies  strongly  with  N.  For  example,  halving  the 
channel  gas  density  at  fixed  pressure  raises  n^/N  by  over  six  orders  of 
magnitude  at  temperatures  T^  <  W^/30.  This  strong  variation  usually  causes 
detracking  condition  (45)  to  be  satisfied  at  channel  temperatures  much 
above  W^/30.  In  hot  air,  for  example,  the  lowest  ionization  potential  is 
Vj  =  9.3  eV  for  NO,  so  that  condition  (45)  is  usually  satisfied  if  the 
channel  temperature  rises  much  above  0.3  eV.  This  temperature  restriction 
limits  the  depth  of  the  channel,  and  thus  limits  the  degree  of  beam-range 
extension  possible. 
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E.  Chemistry  Summary 

The  preceding  analysis  of  gas  chemistry  suggests  that  the  beam  head  is 
likely  to  track  the  channel,  as  long  as  channel  preionization  and  on-axis 
avalanching  are  weak  or  absent.  Although  the  beam  body  may  experience 
detracking  from  higher-order  chemistry,  longitudinal  coupling  causes  the 
body  to  follow  the  head.  The  sign  and  magnitude  of  the  channel  force  is 
thus  most  important  in  the  region  where  coupling  first  dominates. 

As  discussed  in  Sec.  V-B,  we  have  verified  the  change  in  sign  of  the 

tracking  force  at  late  t,  as  a  function  of  various  chemistry  processes. 

19 

More  exact  chemistry  analysis  by  Keeley  shows  qualitatively  similar 
results.  Ve  have  also  performed  dynamical  simulations,  to  be  reported  in  a 
later  paper,  showing  that  the  entire  beam  can  be  pulled  into  the  channel, 
even  when  the  channel  force  is  detracking  in  the  body.  The  first  such 
simulation  was  performed  by  Welch**  using  the  particle  code  IPROP.  The 
dynamical  codes  employ  relatively  simple  chemistry,  and  they  initially 
failed  to  detect  density  tracking  because  was  taken  to  be  constant. 

VI.  STATIC  SIMULATIONS 

A.  Comparison  vith  Analytic  Theory 

We  have  used  two  computer  codes,  VIPER-1  and  SARLAC,  to  compute  the 
average  force  produced  by  a  density  channel  on  a  rigid-rod  beam.  For  this 
purpose,  the  beam  is  not  propagated  in  z  but  is  treated  as  frozen.  In 
dynamical  runs,  which  will  be  described  in  a  later  paper,  the  beam  is 
allowed  to  move,  although  Maxwell's  equations  are  still  solved  using  the 
frozen  approximation,  d/dz  =  0.  VIPER  is  a  linearized  code  limited  to 
small  beam  displacements,  while  SARLAC  is  a  nonlinear  particle  code. 
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To  test  the  nonlinear  theory,  SARLAC  vas  run  using  beam  impact 
ionization  only,  a  density  channel  given  by  Eq.  (14),  and  a  collision 
frequency  given  by 


Vr-e-« 


V 

mo 


1  -  q  ln[(E/N)/(Eo/No>] 


(46a) 


where  E  is  the  total  electric  field  and  q,  v  ,  and  E  /N  are  constants. 

mo  oo 

The  form  of  this  equation  requires  a  normalization  coefficient  given  by 


k(6)  *  1  +  §9-  [ln(8)  +  gl(l)  -  g1(6)]  -  q  In  ,  (46b) 

v^*o  oy 


in  place  of  Eq.  (20d).  Here  E'  equals  1^  divided  by  the  plasma 
conductance.  In  the  results  presented  below,  E'  was  set  for  convenience  to 
the  on-axis  electric  field  strength. 

Simulations  were  performed  using  a  beam  with  a  Bennett  profile  of 
constant  radius  rfa  «  1  cm  and  a  current  Ib  that  rose  linearly  to  10  kA  in 
10  ns  and  then  remained  constant.  The  parameters  8,  p,  yc,  r£,  q,  vmQ,  and 
Eq/Nq  were  varied  over  more  than  an  order  of  magnitude.  In  Fig.  4  we  plot 
the  net  deflection  force,  Ft(0,  as  computed  by  SARLAC  for  several  typical 
runs.  This  force,  which  includes  contributions  from  the  plasma  charges  and 
displacement  currents,  agrees  well  with  the  magnetic  channel  force,  Ftm(C) 
given  by  the  theory,  after  a  few  ns.  Observe  that  the  peak  tracking  force, 
6  Gauss  in  Fig.  4a,  is  roughly  two  orders  of  magnitude  smaller  than  the 
monopole  pinch  force  at  that  point  in  the  beam.  Beam  distortion  should 
therefore  be  negligible.  The  detracking  exhibited  at  early  C  is  an 
electrostatic  artifact  of  SARLAC,  and  is  not  seen  in  VIPER  (or  in  SARLAC  if 
the  beam  head  is  flared.)  Dynamical  runs,  to  be  described  in  a  later 
paper,  show  the  beam  being  pulled  into  the  channel  after  propagating  a 
distance  consistent  with  Eq.  (12). 
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B.  Gas  Cheaistry 

To  examine  the  effects  of  various  chemistry  processes  on  channel 

tracking,  we  performed  a  series  of  runs  using  the  linearized  code  VIPER. 

In  these  runs,  the  beam  current  rose  smoothly  from  zero  to  20  kA,  and  the 

beam  radius  fell  smoothly  from  2.65  cm  to  1  cm,  over  a  time  scale  of  10  ns 

(AC  =  300  cm).  The  radius  variation  approximates  the  flaring  that  develops 

25 

naturally  from  weak  pinching  of  the  beam  head.  A  Bennett  profile  was 
used  for  the  beam,  and  a  Bennett  profile  of  radius  r£  =  3  cm  was  used  for 
the  channel  depth.  The  ambient  gas  density  was  Nq  =  1  atm,  and  the  on-axis 
channel  density  was  N  =0.3  atm. 
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VIPER  employs  an  air  chemistry  model  containing  all  but  the 

convection  term  in  the  continuity  equation  (33).  In  Fig.  5  we  plot  the  net 

channel  force  as  a  function  of  £  as  various  chemistry  processes  are  turned 

on  separately  and  in  total.  Curve  (a)  plots  F  /y  using  beam  ionization 

only;  the  force  is  everywhere  tracking  and  agrees  well  with  the  analytic 

predictions.  Curve  (b)  adds  electron  attachment  to  02  at  a  rate  aflNo  - 

5x10^  s-1;  attachment  weakens  the  density  channel  force  consistent  with  Eq. 

(35),  but  does  not  cause  actual  detracking  until  much  later.  Curve  (c) 

shows  the  effect  of  electron-ion  recombination  with  a  rate  coefficient  &r  - 
-7  3 

10  cm  /s;  recombination  weakens  the  density  channel  force  consistent  with 
Eq.  (36),  but  again  does  not  cause  detracking  until  much  later.  Curve  (d) 
shows  the  effect  of  Spitzer  collisions  which,  for  the  given  parameters, 
degrade  detracking  more  than  attachment  or  recombination.  Curve  (e)  shows 
the  cumulative  effect  of  all  these  processes.  Observe  that  the  channel 
force  remains  tracking  to  (  >  600  cm,  well  into  the  region  of  strong 
magnetic  coupling.  The  entire  beam  is  thus  predicted  to  track  the  channel. 
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Shortening  the  beam  rise  time  increases  the  inductive  electric  field 
Ez-  This  raises  Tg  and  so  that  the  plasma  conductivity  decreases,  until 
the  field  is  so  high  that  avalanching  commences.  For  the  beam  described 
above,  VIPER  shows  avalanche- induced  detracking  if  the  beam  current  rises 
faster  than  30  kA/ns. 

Adding  channel  preionization  in  excess  of  the  overheat  condition  (40) 
produced  strong  detracking  of  the  beam  head  in  the  VIPER  runs.  Much  lover 
levels  of  preionization  produced  weak  channel  forces  that  were  attractive 
or  repulsive,  on  average,  depending  on  £  and  on  the  relative  sizes  of  the 
beam  and  channel.  Preliminary  dynamical  simulations  suggest,  however,  that 
beam  expulsion  occurs  only  if  the  revised  overheat  condition  (45)  is  met. 


VII.  CONCLUSION 


The  analytic  and  numerical  work  presented  here  strengthens  the 
theoretical  base  for  the  density  tracking  force  discovered  by  Welch. ** 
Three  chemistry  properties  underlie  the  effect.  First,  beam  impact 
ionization  is  the  dominant  source  of  ionization,  and  it  produces  a  degree 
of  ionization  that  is  symmetric  about  the  beam.  Second,  the  low-density 
channel  reduces  the  collisional  cooling  rate  of  the  plasma  electrons,  thus 
raising  their  temperature.  And  third,  in  most  gases,  the  elevated 
temperature  raises  the  momentum- transfer  collision  frequency  of  the  plasma 
electrons,  thereby  lovering  the  plasma  conductivity.  The  reduced 
conductivity  lovers  the  plasma  return  current  in  the  channel,  so  that  the 
beam  is  magnetically  attracted  to  the  channel. 

Several  properties  of  density  tracking  make  it  unusually  effective  at 
steering  a  beam.  First,  the  density-channel  force  is  strong,  potentially 
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as  strong  as  the  radial  pinch  force  that  binds  the  beam  together.  Second, 

the  density  force  is  long-lived,  because  it  develops  from  beam  ionization 

and  from  magnetic  monopole  fields,  rather  than  from  magnetic  dipole  or 

electrostatic  fields.  Third,  the  density  force  transitions  smoothly  to  a 

longitudinal  coupling  force  that  causes  the  beam  body  to  follow  the  head, 

regardless  of  the  sign  of  the  channel  force  in  the  body.  Fourth,  the 

channel  force  is  tracking  in  the  beam  head  over  a  vide  parameter  range,  and 

remains  tracking  even  as  the  beam  pinches  to  a  narrow  radius.  This  trait 
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is  especially  important  for  beams  that  are  radius-tailored  to  reduce  the 

grovth  rate  of  the  resistive  hose  instability.  Fifth,  the  density  force 

helps  stabilize  the  hose  instability,  particularly  in  the  beam  head  where 

stabilization  is  most  needed.  And  sixth,  density  tracking  develops 

naturally  out  of  the  need  to  achieve  range  extension.  The  principal 

requirements  for  density  tracking  to  occur  are  low  channel  preionization, 

veak  on-axis  avalanching,  and  a  chemistry  parameter  q  that  is  positive. 

In  short,  density  tracking  should  be  capable,  under  appropriate 

conditions,  of  keeping  a  stable,  relativistic  electron  beam  inside  a 

11-13 

rarefied  gas  channel.  Experiments  designed  to  test  the  theoretical 
predictions  have  verified  the  existence  of  a  robust  density-tracking  force, 
although  they  have  not  yet  demonstrated  range  extension  or  clear  evidence 
of  longitudinal  coupling.  The  latter  issues  require  longer  propagation 
distances  and  more  stable  beams. 
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Appendix:  Longitudinal  Coupling 


Longitudinal  coupling  arises  when  the  beam  and  local  pinch  force  are 
misaligned.  The  misalignment  produces  a  transverse  restoring  force, 
denoted  the  coupling  force,  on  the  beam.  For  small  misalignments,  the 
coupling  force  can  be  computed  by  linearizing  the  pinch  force  about  the 
beam  axis.  For  a  self-similar  beam  in  the  magnetic  regime,  the  coupling 
force  takes  the  form 


F 

c 


eIn  *b  ~  ?n 
rbc  rb  ’ 


(Al) 


where  In  is  the  net  current,  yn  is  its  centroid,  y^  is  the  beam  centroid 
(replacing  yc  in  the  earlier  analysis),  and  e^/r^c  is  the  average  magnetic 
pinch  force.  Here  we  assume  y^>yn  «  where  y^  and  yc  are  measured  with 
respect  to  a  fixed  (density-channel)  axis.  In  addition,  In  is  taken  to  be 
symmetric  about  its  centroid  yn>  and  azimuthal  asymmetries  from  a  channel 
are  treated  separately  as  a  channel  force. 

Ve  previously  found  that  the  magnetic  dipole  fields  responsible  for 
magnetic  coupling  relax  on  a  dipole  decay  length  This  suggests  that 

the  net  current  centroid  centers  about  the  beam  according  to 


(A2) 


Combining  Eqs.  (Al)  and  (A2)  yields  the  approximation 
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eIn 

■  V  rb  «  ■ 


A  dynamical  equation  for  dy^/dz  closes  the  problem. 


14,16 


(A3) 
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* 

To  calculate  the  channel  force,  ve  had  assumed  a  rigid-rod  beam 
parallel  to  the  channel.  Ve  could  then  ignore  the  coupling  force  Fc> 
because  there  was  no  tilt  to  either  the  beam  or  net  current,  3yb/3C  = 

3y  /3C  =  0  *  y.  -  y  .  However,  a  channel  force  that  varies  with  £  soon 
causes  the  beam  to  separate  from  the  net  current  and  to  tilt,  so  that 
yn  *  yb,  3yn/3C  *  0*  and  Fc  *  0.  A  large  beam  tilt,  combined  with  large 
can  produce  a  coupling  force  larger  than  the  channel  force.  In  the 
beam  body,  is  large,  and  coupling  usually  dominates  so  that  the  body 
follows  the  head. 

The  precise  point  at  which  coupling  dominates  a  given  channel  force  is 
complicated  by  the  different  dependencies  of  the  forces  on  the  beam  and 
channel  parameters.  However,  for  a  channel  force  that  transitions  from 
detracking  to  tracking  or  vice  versa,  the  coupling  force  is  likely  to 
dominate  provided  the  transition  occurs  quickly  relative  to  In  the 

case  of  channel  preionization,  condition  (45)  should  be  sufficient  for 
coupling  to  occur. 

An  interesting  observation  is  that  condition  (45)  can  be  met  only  at 
high  beam  currents.  To  show  this,  let  us  use  the  left-hand  side  of  the 
continuity  equation  (33)  to  obtain  the  general  expression 

>  ecnrb2n(Ct)/^siIb(Ct),  (A4) 


depending  on  how  fast  Ib  rises  with  C>  If  ve  now  set  n(Ct)  =  n^,  based  on 
Eq.  (44),  ve  find  that  condition  (45)  can  be  rewritten  as 
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-7  3  —  id  2 

Typical  values  in  air  are  vn  a  10”  cm  /s,  q  =  0.2,  and  s^  a  3x10”  cm  , 

suggesting  that  condition  (45)  is  satisfied  only  for  I^(Ct)  >  25  kA.  If 

the  beam  current  never  rises  to  this  value,  coupling  requirement  (45)  is 

never  satisfied,  and  the  beam  body  need  not  follow  the  head  but  is  free  to 

track  the  channel.  In  principle,  then,  the  body  of  modest-current  beams 

can  track  a  density  channel,  regardless  of  preionization  level. 

In. practice,  hovever,  high  levels  of  preionization  push  Ct  into  regions 
where  higher-order  chemistry  causes  the  density  force  to  become  detracking. 
The  net  channel  force  is  then  everywhere  detracking,  so  that  coupling  is 
irrelevant,  and  the  entire  beam  is  ejected  from  the  channel.  Furthermore, 
high  channel  preionization  will  excite  violent  hose  instability.^ 
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Fig.  2a  —  Density  force  F^ff)  as  a  function  of  channel  offset  yc.  F^  is  given  in  units  of  Gauss  per 
kA  of  plasma  return  current,  calculated  using  Eqs.  (19)-(20)  with  the  following  nominal  parameters: 
rb  =  rc  =  yc  =  1  cm,  q  =  0.2,  6  =  0.1,  and  p  =  1. 
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Fig.  2d  —  Density  force  as  a  function  of  channel  depth  5. 
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Fig.  3  —  The  inverse  of  the  reduced  collision  frequency  ^asa  function  of  electric  field  in  air. 
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Fig.  4a  —  Typical  comparison  plots  of  the  total  numerical  tracking  force  Ft  [solid  line]  and  the 
theoretical  magnetic  force  Fm  [dashed  line]  for  >yc  =  1,  rc  =  2,  q  =  0.2,  b  —  0.1,  and  p  -  1. 
The  lengths  yc  and  rc  are  normalized  to  the  beam  radius,  rb  =  1  cm. 


44 


Ft(G) 


£  (cm) 

Fig.  4b  —  Typical  comparison  plots  of  the  total  numerical  tracking  force  F,  [solid  line]  and  the 
theoretical  magnetic  force  Fta  [dashed  line]  for  ye  =  3,  rc  *  1,  q  -  0.1,  6  -  0.1,  and  p  =  1. 
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Fig.  4c  —  Typical  comparison  plots  of  the  total  numerical  tracking  force,  Ft  [solid  line]  and  the 
theoretical  magnetic  force  F^  [dashed  line]  for  yc  =  0.1,  rc  =  0.5,  q  +  0.4,  5  =  10~5,  and 

p  =  0.2. 
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Fig.  4d  —  Typical  comparison  plots  of  the  total  numerical  tracking  force  F,  [solid  linel  and  the 
theoretical  magnetic  force  F®  [dashed  line]  for  yc  =  0.1,  rc  =  2,  q  =  0.01,  6  =  0.63,  and  p  =  5. 
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